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Abstract 

A progress report on two recent theoretical approaches proposed to understand the 
physics of irreversible fractal aggregates showing up a structural transition from a rather 
dense to a more multibranched growth is presented. In the first approach the transition 
is understood by solving the Poisson equation on a squared lattice. The second approach 
is based on the discretization of the Biharmonic equation. Within these models the tran- 
sition appears when the growth velocity at the fractal surface presents a minimum. The 
effects of the surrounding medium and geometrical constraints for the seed particles are 
considered. By using the optical diffraction method, the structural transition is further 
characterized by a decrease in the fractal dimension for this peculiar class of aggregates. 
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1 Introduction 



Some irreversible fractal aggregates display in nature peculiar structural transitions dur- 
ing their growing processes. Electrochemical deposition experiments and bacterial 
colony growth are two examples of such phenomena. It has been reported that, at a 
certain threshold distance far from the cluster center, these systems exhibit an intrigu- 
ing transition from a rather dense to a more multibranched growth of which little is 
understood. 

Besides the classical diffusion limited aggregation for 'regular' fractals (having fractal 
dimension df ~ 1.7) || §], the Laplacian aggregation model has been proposed which is 
based on the discretized Laplace equation 

4>i,j = \{<t>i+lj + + + > (!) 

with the potential growth site. However, neither of these simple approaches lead to 
explain a complex structural transition during irreversible fractal growth. 

Two alternative attempts have recently been reported in the literature to understand 
the physics behind a structural transition in fractal growth. In the first approach, due 
to Louis et al. || [7], ||, the transition is derived by solving the Poisson equation (on a 
squared lattice) 

V 2 = A 2 , (2) 

which becomes dependent on the potential 4> at two boundaries, the distance L between 
them, and a screening length A. 

The second theoretical approach is due to Wang and Canessa H [UJ and is based on 



the Biharmonic equation in two-dimensional (2D) isotropic defect-free media, namely 
V 2 (V 2 w) = . (3) 
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By discretizing this equation, these authors have shown that a structural transition can 
also be a consequence of the different coupling of displacements u within the pattern 
formation. 

Both of these approaches describe a similar class of complex topology from different 
perspectives and on different systems. However, the Poisson and Biharmonic growth 
models identically relate the growth probability for each site to the local field as 

P « <p or V 2 m , (4) 

These assumptions phenomenologically follow the Laplacian (or dielectric breakdown) 
model. 

It is the aim of this progress report to review the Poisson and Biharmonic approaches. 
We shall focus on these two model for irreversible fractal aggregates showing up a struc- 
tural transition and shall consider the effects of the surrounding medium and geometrical 



constraints for the seed particles. Following recent ideas by Perez-Rodriguez et al. |TT 
we shall also shown how the optical diffraction method leads to characterize the struc- 
tural transition by a decrease in the fractal dimension for this class of aggregates. 

2 The Poisson Growth Model 



The origin of screening in the Poisson model lies in the presence of free charges and 
leads to a rich variety of patterns ||. As compared to the Laplacian model of Eq.([j]), 
this model introduces a new length scale, i.e., A, and a nontrivial dependence on the 
boundary conditions that are responsible for a structural transition on growing. 

Fig. 1(a) shows an example of the Poisson model in planar geometry by adding 4051 
particles to the cluster and using (j) 1 = 5 x 10~ 7 , <p° = 1 and A -1 = 10. From this figure 
it can be seen that the Poisson patterns obtained can have a fractal character (at scales 
shorter than A), to then display a structural transition. 



Louis at. H have analytically demostrated that this structural transition can be 
characterized by a change in the sign of the electrostatic field at the surface of the 
aggregate besides the minimum of the growth velocity (to be shown later). However, 
Wan Wei || has shown that such a transition may be altered by the existence of a 
critical field. 

It is important to note that within the Poisson model several fractal branch might 
also grow when attaching more than one particle at each computer step. For details see 
Refs.ll,!. 

3 The Biharmonic Growth Model 

Fig. 1(b) shows the final stage of a Biharmonic fractal displaying features (in circular 
geometry) of a transition. Below the transition point r^, the fractal dimension approaches 
the value for Laplacian growth within error bars. To generate this class of patterns one 
solves numerically the Biharmonic Eq.(|3|). Its discrete form on the lattice site, 

yields the following expression 

Ui-2,3 + 2ui-!j-i - Sui-ij + 2ui-i }j+1 + Uij-2 - Suij-i + 20u itj (5) 
—8uij + i + + 2ui+i,j-i — 8ui + ij + 2'Uj + ij + i + Ui + 2,j — 



requiring the values of the normal derivative for the order parameter u |fL3| . For the 
sake of simplicity one sets the derivative boundary condition (necessary along the radial 
direction) equal to zero. 

Throughout calculations one uses lattice sites enclosed within a circle of (normal- 



ized) radius r = \Ji 2 + j 2 such that u° and u l are unity and zero at the outer circular 
boundary and the inner growing Biharmonic aggregate, respectively. The seed particle 
can be placed either centered or distributed under different geometrical constraints as 
discussed later. The procedure for growing fractals then follows standard techniques 



H, till solutions of the discretization of Eq.(Q) converge. Aggregates then stochastically 
grow under the relation between the growth probability P (at the grid site and u 

as in Eq.(|). 



4 Growth Velocity 

Results for the grow velocity v along the y-direction for the Biharmonic model are plotted 
in Fig. 2. For comparison, also included in this figure are the results for Poisson growth 
assuming v to be proportional to the field | fcj — (p 1 | by following Ref.||. 

From Fig. 2, it can be seen that the structural transition coincides with the fact that 
v on the growing surface presents a minimum. The same is true for Biharmonic growth 
(even independently of how one relates the probability Pij to u(i,j) 0). For Laplacian 
growth this phenomenon does not appear because the trend is to generate a single tip 
at faster velocity than in the cases of Poisson or Biharmonic growth. 

Biharmonic patterns below re, (i.e., within the dense region of Fig. 1(b)), does not 
become that dense as for Poisson of Fig. 1(a). This effect can be understood from the 
velocity plot since for Poisson growth the transition occurs at smaller velocities than for 
Biharmonic growth. Henceforth, a Eden-like pattern can be generated due to screening. 
Both Poisson and Biharmonic models lead to growth velocities along the ?/-direction 
with parallel slopes. 

5 Effects of Surrounding Medium and Geometrical 
Constraints 

Within the Biharmonic model, the effects on fractal growth of the surrounding medium 
and of different geometrical constraints (or boundary conditions) for the seed particles 
has been analysed in Ref . [1(J by following the rj (or dielectric breakdown) model || [14| 



5 



and assuming an extention of Eq. (|4|) as 

p. . = 1 V 1 (Q) 



where the sum runs over nearest neighbor sites and rj > 0. 

Figure 3(a) shows typical numerical results by using n = 2 for a cluster of about 
1000 particles. In the limit rj — > oo, a 'transition from slow to faster growth' is found in 
such a way that one end of the needle-like structure presents greater growth probability 
than the other. Such structures become dendritic below and above a transition point. 
On decreasing the value of n, one finds a transition from dendritic-to- compact growth 
such that the inner region of the aggregates becomes denser. If rj — ► 0, the growth 
probability becomes purely random and independent of the Biharmonic field -as in a 
Eden-like pattern. 

The effects of geometrical constraints for the seed particles has also been studied in 
Ref. |10] by fixing the value of rj and using two different geometries for the seed particles, 



namely circular and linear. Figure 3(b) shows a Biharmonic fractal under the constrains 
of a circular (empty) area of seed particles. As can be seen a structural transition still 
survives independently of such a simple configuration adopted. 

6 Optical Diffraction Analysis 

The effect of the structural transition on the diffracted intensity and, consequently, on 
df follows by considering the fractal to be composed of N identical and similarly oriented 



particles on the plane x — y fL5| , jig] . The position of their centers of mass is given by 
— (x n , Un)-, where n — 1, N. In this system, the Fraunhofer diffraction pattern for 
the fractal is derived by assuming that each particle corresponds to one aperture. 

The form factor, corresponding to the intensity scattered by one 'aperture', is deter- 
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mined by the integral of the diffraction amplitude. The structure factor is 
1 N 

S(k)=\-Y,e-^\ 2 , (7) 

n=l 

where k is the component, parallel to the xy-plane, of the scattered wave vector whose 
modulus is k — ^sin# « ^-9, with 9 the (small) angle that the scattered wave vector 
makes with the z-axis. A is the wavelength of the incident light. If ka < 1 (a being 
the size of an elementary particle), the form factor is practically unity and the light 
distribution in the diffraction pattern is given by the structure factor such that the 
normalized diffraction intensity is 

J(k)«S(k) . (8) 

For jr < ka < 1 (L being the size of the whole aggregate) the expected value of the 
intensity -or, alternatively, of the structure factor S(k)- is 

< J(k) >= j d 2 Re' lkK { J d 2 R p(R )p(R + R ) ~ iT") . (9) 

This relation for fractal aggregates obeys a power law behaviour || , where the exponent 
a is related to the Hausdorff dimension df = d — a with d the Euclidean space dimension. 
This variation, in turn, leads to the power law behavior of the diffracted intensity as a 
function of the wave vector 

/(k) ~ kr d * . (10) 

Information about df can be obtained at low k values (such that j- <^c| k x a |<C 1) 
from the averaged intensity defined by the expression 

< I(k) >= — <f d(f)I(\i) ; k x = kcoscj), k y = ksmcf), (11) 

As seen in Fig.4, this quantity for Poisson fractals has a power law behavior as a function 
of the modulus k. However, the structural transition during fractal growth leads to a 



change in df. The slopes of this plot represent a decrease of the fractal dimension for 
Poisson aggregates. Similar results have been found in Biharmonic fractals [11 . 



7 Discussion 

In both models studied the global influence of a growing pattern to the growth probability 
for each lattice site under a certain power-law form is set phenomenologically as in Eq. ([|) 
or (P). However an important difference to note is that within the Biharmonic model, 
iterative procedures are carried out around thirteen next nearest neighbors and not on 
four as within the Poisson and Laplacian models. Thus the formation of connected 
patterns within the Biharmonic equation becomes non trivial and more involved. 

Within the Poisson model, screening lies in the presence of free charges. The physical 
relevance for the Biharmonic equation might follow from elasticity theory (the deflection 
of a thin plate subjected to uniform loading over its surface with fixed edges), the steady 
slow 2D motion of a viscous fluid or the vibration modes in the acoustic of drums. Besides 
this, a higher order differential equation containing the Biharmonic term also appears 
in the study of kinetic growth with surface relaxation (see, e.g., recent work in [|T^] and 
references therein). 

The influence of ramified Biharmonic patterns on the growth probability for each 
lattice site is on the type of pattern obtained and not on locating the transition. Sev- 
eral branches may also develop within the Poisson and Biharmonic growth models by 
attaching -simultaneously, and stochastically- more than one particle at each time step. 
For Poisson growth the structural transition is the result of including many-body con- 
tributions via screening in a sort of mean field approach. In the Biharmonic model, long 
range coupling appears naturally as a consequence of discretizing Eq.(0) (as given in 
Eq.(B])). For planar symmetry, the transition point within both models appears when 
the growth velocities exhibit a minimum. This implies that Poisson and Biharmonic 
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growth describe a similar class of complex structural transition phenonema from two 
different perspectives and on different systems. 

The structural transition point can, to a good approximation, be estimated from the 
continuous limit of Eq.(|3|) in cylindrical coordinates. Canessa and Wang [TO] have shown 



that the transition point depends on the system size L only and occurs approximately 
at a distance (in lattice units) about 60% far from the seed particle indicated by the 
circle arrows in Fig. 1(b). This prediction is in accord with their numerical simulations 
for 77 = 1. 

By tuning 77 — » in Eq.@, the structural transition corresponds to a 'dense-to- 
multibranched transition' whereas for rj — > 00 one founds a 'transition from slow to 
faster growth'. On decreasing 77 from infinity to zero, Biharmonic fractals become denser 
as for Laplacian growth. For large 77, the needle-like structures grow faster at one end 
and, in the overall, remain denditric. Furthermore the transition point is independent 
of geometrical constrains adopted for the seed particles. 

Results of optical diffraction have enabled to identify and relate changes in df of ag- 
gregates to variations in the diffracted intensity as a function of the wave vector. Such 
findings may be experimentally confirmed. In fact, the averaged intensity < I(k) > 
might be determined by using the set up as the one described in [I7|. Therein, a photo- 



multiplier connected to a multichannel analyzer records I(k x , k y ) and the displacement 
of the photomultiplier is controlled by a high-precision motorized micrometer. Then, 
after scanning the diffraction patterns displaying a structural transition, the average of 
intensity over concentric circles (< I(k) >) might be obtained. 
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ure captions 



Fig.l: Different types of irreversible aggregates displaying the transition, (a) 
Poisson fractal in planar geometry, (b) Biharmonic fractal in circular geometry. 

Fig. 2: Growth velocity v along one direction for two different Poisson fractals 
(grouped as I) and for three different Biharmonic fractals (grouped as II). 

Fig. 3: (a) Biharmonic fractal in circular geometry using in Eq.(^) rj — 2 to mimic 
surrounding medium effects, (b) Biharmonic fractal under the constrains of a cir- 
cular (empty) area of seed particles. 

Fig. 4: Log-Log plot of the angle-averaged < I(k) > against k for a typical Poisson 
fractal. Full and open squares, represent results for this class of cluster before and 
after the structural transition. 
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